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The combined effect of both nonmagnetic and magnetic impurities on the superconducting transition 
temperature is studied theoretically within the BCS model. An expression for the critical temperature 
as a function of potential and spin-flip scattering rates is derived for a two-dimensional superconductor 
with arbitrary in-plane anisotropy of the superconducting order parameter, ranging from isotropic s-wave 

CI 

to (i-wave (or any pairing state with nonzero angular momentum) and including anisotropic s-wave and 
mixed (d + s)-wave as particular cases. This expression generalizes the well-known Abrikosov-Gor'kov 
formula for the critical temperature of impure superconductors. The effect of defects and impurities in 
O j high temperature superconductors is discussed. 

2 ■ 1. INTRODUCTION 

The mechanism of high-T c superconductivity still remains unknown. It is generally believed that elucidation of the 
symmetry of the superconducting order parameter A(p) in high-Tc superconductors (HTSCs) could narrow the list 
of pairing mechanisms debated. Since the CuC>2 layers are thought to be responsible for the superconducting pairing 
* ■ in HTSCs, the in-plane symmetry of A(p) is of primary interest. However, in spite of strong evidence provided for 
d x 2_ y 2 in-plane symmetry of A(p), experimental data are somewhat controversial |[-16| (though the recent research 
seems to resolve the contradiction in favor of a d 3 .2_j / 2-wave |T^]). The situation is complicated by orthorhombic 
distortion of CuC>2 layers in some HTSCs, resulting in probable admixture of s-wave component to otherwise pure 
d-wave A(p) |l^JT^]. Since HTSCs differ structurally, the form of A(p) may appear to be material-dependent. 
I ' One indirect way to distinguish a pure d-wave from a highly anisotropic s-wave or a mixed (d + s)-wave is to 
. study the response of HTSCs to intentionally incorporated impurities or radiation- induced defects. Depending on the 
symmetry of A(p), clear differences are predicted for the variation of experimentally accessible characteristics such 
as the critical temperature T c , the superfluid density, etc pp|-|3^]. For example, the rate of T c degradation by defects 
and impurities in a two-dimensional superconductor is determined by the value of the Fermi surface (FS) average 
(A(p)) ps [§0| and should be different in d-wave superconductors with (A(p)) fs = and anisotropic s-wave or mixed 
(d + s)-wave ones with (A(p))i?s 7^ 0, the specific value of (A(p))ps being dictated by the degree of A(p) anisotropy 
OS . or by the relative contributions of d-wave and s-wave components to A(p). 

Numerous experimental studies give evidence for T c degradation by impurity doping |]33| -^8| or radiation damage 
j)|-[58f of HTSCs, the two effects being remarkably analogous if the dependence of T c on the in-plane residual resistivity 
Po is considered f39f . Impurity-induced scattering of charge carriers in doped HTSCs and their scattering by displaced 
1 host atoms in irradiated HTSCs are believed to be the main reasons for the suppression of superconductivity and the 
increase in Po ||,|§,||,||J5|J5|,||,|| . 

The comparison of experimental curves T c (po) with theoretical ones reveals that the observed reduction of 
T c by impurities and radiation defects is more gradual than predicted theoretically for d-wave superconductors 
~l|,^0 53 5qj57| , |59| - |6]]] . A critical value of at which T c = ranges from 200 pVLcm to 1500 /ificm depending 
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on the type of disorder and the kind of HTSC material [ Pj4T|^,P ,^|,p|,p|j55|,|57| H , while for a d-wave supercon- 
ductor with T c w 100 K the theory gives ~ 50 /iOcm ]3l|j59| -|6l|" To reconcile the experimental findings with the 
d-wave symmetry of A(p) in HTSCs, a number of suggestions have been made, including the anisotropy of impurity 
scattering in the momentum space J62| , |63[ , an " intermediate" (between Cooper pairs and local bosons) state of paired 
electrons |54j , a depletion of the hole density due to the oxygen vacancies in the Cu02 planes [ p5[, an anomalously 
small value of the plasma frequency J56|;Q , the spatial variation of the order parameter J66[ , etc ]67| . 

Another way is to abandon the d-wave hypothesis in favor of anisotropic s-wave or mixed [d + s)-wave models 



[|27j,[28 J30| 3lj,|68|] (it must be emphasized that if the relative weight of the isotropic s-wave component in the mixed 
(d + s)-wave A(p) is large, then the symmetry of such an order parameter can in fact be viewed as the anisotropic 
s-wave one j53fl). However, while the initial slope of experimentally observed T c (po) curve in HTSCs can actually 



be explained by anisotropic s-wave symmetry of A(p) [pl| , p3U 57|, the theory faces problems when explaining the 
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complete suppression of superconductivity at a finite value of p . Indeed, the theory predicts that T c of 

a non-<i-wave anisotropic two-dimensional superconductor doesn't vanish at a certain critical value of p% (as it does in 
the case of pure d-wave symmetry of A(p)), but instead asymptotically goes to zero as po increases. This contradicts 
the experiments mentioned above. Besides, the experimentally observed form of the T c (/?o) curve is usually close to 
linear ]42] , ^ , p7| , p r 2| , p 7 3|j5^j57| |59| , while the theory predicts a positive curvature of the T c (po) curve in a non-d-wave 
superconductor . 

Note, however, that experimentally determined values of po reflect the contribution from different scattering chan- 
nels, while a theoretical analysis of T c degradation by defects and impurities in HTSCs is usually made for the specific 
case of spin- independent scattering potential 20| 2^, |25| , p7| |2 53 5^, |57| , |60 62,63-167]]. Meanwhile a lot of experiments 
give evidence for the presence of magnetic scatterers (along with nonmagnetic ones) in non-stoichiometric HTSCs, e.g., 



For example, the oxygen vacancies or excess oxygen 



in oxygen-deficient, doped or irradiated samples [ |36| , |45[ |6 

atoms carry or induce local magnetic moments and hence play a role of paramagnetic centers 
doping by Zn induces local magnetic moments residing probably on the nearest-neighbor Cu sites 




1|. Furthermore, 
frail. This is 



supported by studies of nonmagnetic impurities in Heisenberg antiferromagnets 



and by numerical calculations 



within the two-dimensional t — J model |82| . Besides, there are intrinsic (host) magnetic atoms in some stoichiometric 
HTSCs, e.g., in GdBa 2 Cu 3 7 !§. 

Though the estimated moment-carrier exchange energy J may appear to be too small to solely account for suppres- 
sion of superconductivity in disordered HTSCs [Q , an intriguing possibility of the magnetic pair-breaking scattering 
as the common origin of the significant decrease in T c remains [§5|. Moreover, it is suggested [§6| that interaction be- 
tween Cooper pairs and localized magnetic moments in " optimally" doped HTSCs leads to a depression of T c relative 
to its "intrinsic" value. Hence, among other things, the understanding of the role of magnetic scattering in HTSCs is 
important from the viewpoint of search for new materials with higher T c . 

Since, first, there exist two channels of carrier scattering by magnetic impurities (potential and spin-flip ones) 
and, second, in general both magnetic and nonmagnetic scatterers are present in HTSCs, there is a need for a 
theoretical model which could describe the effects of nonmagnetic and magnetic scattering on equal footing. For an 
isotropic s-wave superconductor this is the Abrikosov-Gor'kov theory [§7| which predicts a rapid T c suppression by 

in accordance with the Anderson theorem |8^] 
IQI seems to be resolved in favor of the standard 



magnetic impurities and inscnsitivity of T c to nonmagnetic scattering 
(a discussion about the validity of the Abrikosov-Gor'kov approach |8{| 
Green's functions technique). However, the Abrikosov-Gor'kov formula for T c versus scattering rate is not applicable 
to anisotropic superconductors, no matter what the specific symmetry of A(p) is (d-wave, (d + s)-wave, anisotropic 
s-wave or somewhat else). Hence, this formula cannot be used to describe impurity effects in HTSCs. On the 
other hand, theoretical considerations of impurity scattering in anisotropic superconductors are commonly restricted 
to nonmagnetic scatterers only |20|pl| , |23] |25|,^] [32 60 6^j64| . Such a status of the theory of impure superconductors 
results in situations when the experiments on the T c reduction by impurities or radiation-induced defects in HTSCs are 
compared with either the Abrikosov-Gor'kov formula for T c of isotropic s-wave superconductor containing magnetic 



impuritie_s_ 
only] 41 
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J ,8 



|r?1,5ll , |70| or the formula for T c of anisotropic superconductor but containing nonmagnetic impurities 



52-55]. In the latter case an a priori suggestion is often made about pure d-wave symmetry of A(p) 
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In a recent study [|9l| , the Abrikosov-Gor'kov theory has been generalized to the case of a multiband superconductor 
with arbitrary anisotropy of interband order parameter and arbitrary strength of magnetic and/or nonmagnetic 
impurity scattering. Note, however, that the concept of multiband superconductivity (arising, e.g., from Cu02 
plains and Cu-0 chains) is hardly probable to be applicable to HTSCs. Although it was argued in |)l] that the 
mathematical formalism had been proven to be the same for a multiband superconductor and a superconductor with 
a general angular anisotropy of the order parameter ]p2[ , an explicit formula for the critical temperature of a one-band 
anisotropic superconductor containing both nonmagnetic and magnetic impurities was not offered in |9l| . 

The goal of this paper is to work out a theoretical framework for a description of combined effect of nonmagnetic 
and magnetic scatterers on T c of a two-dimensional superconductor with anisotropic A(p) (preliminary results have 
been presented in J9||). We seek to obtain a rather simple (free from needless theoretical complications) Abrikosov- 
Gor'kov- like formula for T c which included physically meaningful parameters and could be compared with available 
experimental data. Within the weak coupling limit of the BCS model and without specifying the microscopic mech- 
anism of superconducting pairing, we derive the expression that relates T c to relaxation rates of charge carriers by 
nonmagnetic and magnetic scatterers, as well as to the numerical coefficient \ — 1 ~~ (A(p))^s/(A 2 (p))^ s which is 
a measure of the degree of in-plane anisotropy of A(p) on the FS. The range < x < 1 covers the cases of isotropic 
s-wave (x = 0), d-wave (x = 1), anisotropic s-wave (0 < x < 1)? an d mixed (d + s)-wave (0 < x < 1) symmetries of 
A(p). In two particular cases of (i) both nonmagnetic and magnetic scattering in an isotropic s-wave superconductor 
(x = 0) and (ii) nonmagnetic scattering only in a superconductor with arbitrary anisotropy of A(p) (0 < x < l)i our 
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expression for T c reduces to the well-known formulae [ |20| , p7| . 

The paper is organized as follows. The BCS model for an impure anisotropic superconductor containing both 
nonmagnetic and magnetic scatterers is described in Sec. 2 along with the theoretical formalism. The expression for 
the critical temperature as a function of potential and spin-flip relaxation times of charge carriers is derived in Sec. 3 
for a superconductor with an arbitrary degree of the order parameter anisotropy. The results obtained are discussed 
in Sec. 4. In Sec. 5 concluding remarks are given. 



2. MODEL AND FORMALISM 

Within the framework of the BCS model, the Hamiltonian of a superconductor containing both nonmagnetic and 
magnetic impurities is as follows 

H = ^£(p)a+ r a p(T + J2 U(p,tT;p',a f )a+ r a ptul +Y, V MKAl"-i>'l«p'V 0) 

p.u p.p',a,(r' P,p' 

where the operator a+ CT (a p(T ) creates (annihilates) an electron with the quasimomentum p and the spin projection on 
z-axis cr = f or j, £(p) = e(p) — fx is the (spin independent) quasiparticle energy measured from the chemical potential 
/i, U(p, a; p', cr') is the matrix element for electron scattering by randomly distributed impurities (defects) from the 
state (p',cr') to the state (p, cr), and T^(p,p') is the BCS pair potential. 
Let the following sum 

U{r) = U n (r)+U m (r) (2) 

be the total interaction between a conduction electron at a point r and all impurities present in the sample, U n (r) 
and U m (r) being the interaction components due to nonmagnetic and magnetic impurities respectively: 

U n (r) = ^U„(r - R a ) , U m (r) = y^u m (r - R^) , (3) 

a (3 

where u n {r — R Q ) is the interaction between an electron at r and a nonmagnetic impurity at R Q , while u m (r — R/3) is 
the interaction between an electron at r and a magnetic impurity at R^. Since, in general, magnetic impurities give 
rise to both potential and exchange scattering [^tJ, one has 

u m (r - Rp) = <*(r - R^) + u«(r - R^), (4) 

where w„*(r — R/3) is the spin-independent potential component, and 

u»(r-R /9 ) = J(r-R /3 )sS / , (5) 

is the exchange interaction. Here J(r — R^) is the exchange energy, Sp is the spin of magnetic impurity located at 
R/3 and s = a/2 is the operator of electron spin (the three components of a are the Pauli matrices a x , <r y , a z ). We 
shall assume that orientations of the paramagnetic spins S/3 are fixed and remain unchanged upon electron scattering 
(taking into account dynamic transitions of the impurity spin between 2S + 1 magnetic sublevels has a little effect on 
the results @). 

The electron interactions u n (r — R a ) and ii^'fr — R/3) with nonmagnetic and magnetic impurities respectively are 
spin-independent and hence contribute to the scattering matrix elements U(p, cr; p', u) only. In its turn, exchange 
interactions of electron with magnetic impurities, u^(r — R/3), result in both spin-conserving and spin-flip scattering 
events and hence contribute to U(jp, a; p', a) as well as to U(p, cr; p', — cr). Let us write down the matrix element 
U(p, cr; p', cr') as 



U(p,a;p',a') = Ui(p,p',o-)6 a <T> + U 2 {p, p', er)cv ,-„> , (6) 



where 



U^p'^a) = Wn (p,p')$>- 4(p - p ' )R « + ^(p^O^e-^P-P')^ + ij(p,p07aE e " i(P_P ' )R " 5 l ' 

/3 S3 

U 2 (p,p',a) = ij(p,p')^ e - j(p - p ' )R ^(^- ? 7^) ■ (7) 
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Here u n (p, p'), u^ >t (p, p'), and J(p, p') are the components of the matrix element for electron scattering by an isolated 
impurity; 7 CT = +1 and -1 for a — f and J. respectively. We do not consider the direct effect of nonmagnetic disorder 
on the magnetic pair breaking (this effect has been studied in Refs. @J9^] for the case of an isotropic s-wave pairing). 

In order to account for anisotropy of the superconducting state, we assume a factorizable phenomenological pairing 
interaction V(p,p') of the form (see, e. g., p0[|) 

F(p,p') = -W(n)0(n'), (8) 
where Vo is the pairing energy, n = p/p is a unit vector along the momentum. Then the order parameter A(p) is ]2C| ] 

A( P ) = -£>(P>P')(fi- P 'iVT> = Ao0(n) ' (9) 
p' 

where Ao depends on the temperature. The function 4>(n) specifies the anisotropy of A(p) in the momentum space 
(e.g., </>(n) = 1 for isotropic s-wave pairing). We assume that </>(n) is temperature independent. 

The self-consistent equation for A(p) can be derived by means of Green's functions technique (see, e.g., fl96fl ). We 
define the normal and anomalous temperature Green's functions 



G(p,cr;p 
F(p,cr;p 
G(p,cr;p 
F(p,cr;p 



*V) = -(^a pCT (r)a+ v ,(0)), 

Sir) = (TVal p _ CT (r)a+ CT ,(0)), 

a';r) = - (T T d+ p _ a (r)a_ p ,_ a ,(0)) , 

ct'-t) = (T T a p Jr)a_ pl _ a ,(0)), (10) 



and their Fourier transforms G(p, cr; p', cr'; w), F(p, <r;p', a'; u), G(p,a;p',a';ui), F(p, a; p', a'; ui), where angular 
brackets stand for the statistical averaging with the Hamiltonian ([!]) , the symbol T T denotes the time ordering, r is 
the imaginary time, and to — ixT(2n + I) are Matsubara frequencies (we set h = kg = 1 throughout the paper). 
It is convenient to introduce the matrix Green function G(p, cr; p', a'; lj) in the Nambu representation: 



1 --F(p,cr;p',cr';w) G(p,cr;p',cr';w ' 



We stress that G(p, cr; p', a'; oj) is nondiagonal in spin space (since there is spin-flip scattering of electrons on magnetic 
impurities) as well as in momentum space (until averaged over impurities coordinates). The matrix equation for 
G(p, cr; p', <t'; cu) can be written as 

Go X (p, a; k, A; w)G(k, A; p', a'; lo) - U(p, a; k, A)G(k, A; p', cr'; u) = i<5 p , p ,<W, ( 12 ) 

where Go(p, cr; k, A; w)=Go(p, cr, w)<5 P! kC'(T,A is the Green function of a clean sample, 

G ( P ,cr,c)-- w2 + e(p)2 + |A(p)|2 7ffA , (p) iw _ e(p) J> (13) 

I is the unit matrix 2x2, and the matrix U(p, cr; k, A) that describes the effect of impurity scattering has the form 



The summation over repeated indices in Eq. ( |12[ ) and below is implied. We note that Go(p, cr, uS) depends on a 
through 7o-. 

In order to avoid needless mathematical complications and to express the final results in terms of as few parameters 
as possible, we make several simplifying assumptions: (i) we consider the short-range scattering potentials, so that 
the matrix elements w„(p, p'), u£?*(p>p')> ^(P>P') are momentum-independent and equal to u n , u^ 4 , J respectively 
(s-wave impurity scattering); (ii) we treat the impurity scattering in the Born limit; (iii) we restrict the momenta of 
the electron self-energy and BCS pair potential to the FS. 

After averaging the Eq. ([12j) over impurity configurations and directions of impurity spins one has 
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where 



G(p, a; p , a ; w) ) = G{p,a,uj)S p ^S a ^, (15) 

imp 



&- l (p,*,w) = G-\p,a,w)-M(p,a,uj), (16) 



M(p,a,w) = (f>(p,a;k, A)G(k,A,^)?7(k,A;p,a; 

(c„| W „| 2 +c m |un 2 +c m |<-| 2 )EkG(k,a, W ) (c n K| 2 + c m K*| 2 -c m |^| 2 )E k ^(k,^w) 

(c>J 2 +C m |< f -C m |<?| 2 )£ k F(k,CT, W ) (Cnl«n| 2 + CmK*| 2 + Cml«"| 2 )EkG( k . ff . w ) 



(17) 



Here c„ and c m are the concentrations of nonmagnetic and magnetic impurities respectively, and we have designated 
\ u m\ 2 — \ J\ 2 S(S + l)/4. Note that | 2 includes contributions from both spin-flip and spin-conserving scattering 
of electrons due to their exchange interaction with magnetic impurities, Eq. (^). This is because the matrix element 
of the spin-conserving exchange scattering depends on the orientation of electron spin through j a , see Eq. (Q). The 
coefficient c m \u e ^\ 2 in Eq. ([Tt]) enters into the factors at the normal and anomalous Green functions with opposite 
signs, while the coefficients c n |w n | 2 and c m |u^'| 2 , which are due to the scattering by nonmagnetic impurities and the 
potential component of the scattering by magnetic impurities respectively, appear in Eq. ( p"7| ) with the same signs. 
Making use of Eqs. (|13|), (|l|), and (jl7|), we obtain 



G(p,a,u,A) = G ( P) <T,uj',A u ), (18) 

where 



-■' = -• - '■«>.. 1 - '■-I'd 2 '-<>;,; 2 )E ^ 2+$ T (k j| ( Jl^ k) | 2 ■ < 



a w ( P ) = A( P ) + ( C „| M „| 2 + c m |<f - Cm io 2 ) E ^+^jt) + |A^(k)p ■ (20) 

Then one has from Eq. (^|): 

A (P)= -T£i:np.rt „ aH ,,^;, (p , )r 

For further considerations it is convenient to express the coefficients c n |M n | 2 , c m |M^*| 2 , and c m |it^| 2 in terms of 
electron relaxation times r„, t^ ', and r^f for scattering by nonmagnetic impurities, potential scattering by magnetic 
impurities, and exchange scattering by magnetic impurities respectively: 

- = 27rc„K| 2 7V(0), _L = 2 7 rc m |< t | 2 7V(0), -L = 2nc m \u%\ 2 N(0), (22) 

where A^(0) is the density of electron states at the Fermi level. The electron relaxation time r m due to magnetic 
impurities is given by the expression 



. - vot + T ex ' ( 23 ) 
m Tm 1 m 



while the total electron relaxation time r due to all impurities present in the sample can be found as 

11111 1 

- = - + - = - + ~FTt+^- (24) 

' 'n 'rn hi T m ' rn 

We note that Eqs. (^) allow one to express the final results in terms of three relaxation times (r„, r^f*, r^f ) instead 
of a large number of unknown parameters such as impurity concentrations and scattering matrix elements. Besides, 
the relaxation times are associated with the residual resistivity. This facilitates a comparison between the theory and 
experiment. 
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3. CRITICAL TEMPERATURE OF IMPURE ANISOTROPIC SUPERCONDUCTOR 



The critical temperature T c can be found from Eqs. (p_9[), (|2C|), and ( |2lD as the temperature at which the order 
parameter goes to zero, i.e., Ao — > in Eq. (0j. Setting |At,;(k)p = in the denominators of Eqs. ( |i9| ) and ( pOf ) and 
taking Eqs. ( p2| ) into account, we have at T — > T c : 

w' = w + ~ (l/r„ + l/r,r + l/<f) « S n(w), (25) 



A^p) = A(p) + ((l/r„ + l/r r r - l/r«) (A^p)}^, 



where the angular brackets (~.)fs stand for a FS average: 



(...) 



FS 



FS y ' |»£(P)/0P| / 7f5 I^(P)/5 P 



(26) 



(27) 



Substituting Eqs. ( p5| ) and ( p6| ) in Eq. (|2~l|), setting |A w (p')| 2 = in the denominator of Eq. (|2l|), and taking Eqs. 
(0) and (0) into account, we have 



1 = T I 

A c 4- M + i(i/ Tn + i 



M + !(l/r n + l/7£* + l/T-) 



(0 2 (n)) FS + (0(n)) 



FS" 



l/ r n + l/rff* - 
2 (M + 1/r-) 



(28) 



where A = VoN(0) is the electron-boson coupling constant. The equation for the critical temperature T c q in the 
absence of impurities (i.e., at l/r n = I/t^ ' = = 0) reads 



\ =^r cO (0 2 (n)) F s^^| 



A 



(29) 



Following the standard procedure, we obtain from Eqs. (|2S]) and ( |2 9| ) the equation for the critical temperature T c as 
T c0 \ _ ^ 1 



In 



^E 



H + 1(1/^ + 1/7^+1/7^) 



J_ /, /T , ! / , , /T e^ _ Wri&s 1/rn + l/rr-l/r^ 
2M (l/r„ + l/r m +l/r„J 2 (M + 1/t~) 

(30) 



At this stage it is convenient to introduce the coefficient x 01 anisotropy of the order parameter on the FS ]20|,pl 



<^(„)) 



= 1 



FS 



(A( P )) FS 
<A»(p)> 



FS 



(31) 



For isotropic s-wave pairing we have A(p) = const on the FS; therefore, (A(p)) 2 ?s . = (A 2 (p))ps, and x = 0. For 
a two-dimensional superconductor with d-wave pairing we have x = 1 since (A(p))jr<j = 0. The range < x < 1 
corresponds to anisotropic s-wave or mixed (d + s)-wave in-plane pairing. The higher the in-plane anisotropy of A(p) 
(e.g., the greater the partial weight of a d-wave in the case of mixed pairing), the closer to unity is the value of X- 

Note that x = 1 holds not only for d-wave pairing state, but also for any pairing state with angular momentum 
I > 0, e.g. for p-wave state (I = 1), see Eq. ((3^). In its turn, the range < x < 1 generally corresponds to mixing of 
s-wave state with some higher angular harmonic state. Hence, while this paper focuses primarily on s-wave, d-wave, 
and (d + s)-wave states, one should keep in mind that the results obtained are more general and may be applied to 
superconductors with other symmetries of the order parameter. 

Making use of the definition ( |3l|) and the formula |)7j 



E 

fe=0 



= - #(sb), 



where ^Sf is the digamma function, we obtain from Eq. (|3C|): 



.n(f| = n-w 



1 



-5 



2 4irT r 



1 

pot 



-* ( 2 



(32) 



(33) 
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In two particular cases of (i) both nonmagnetic and magnetic scattering in an isotropic s-wave superconductor (x = 0) 
and (ii) nonmagnetic scattering only in a superconductor with arbitrary in-plane anisotropy of A(p) (1/r^f = 1 /t*J, = 
0, < x < 1), the Eq. (R3h reduces to well-known expressions |8^j20|| 



In 



2ttT c t- 



(34) 



and 



In 



~T~e 



1 



4ttT c t„ 



(35) 



respectively. 

Now let us consider the limiting cases of low and high impurity concentration (T c0 — T c « T c0 and T c — > 
respectively). At 1/47tT c ot„ << 1, 1 / 'AnT^r^ << 1 and 1/AttT c qt^ << 1 (low impurity concentration) one has 
fromEq. ©: 



T, 



rO 



m 71" 



1 



1 

pot 



1-X/2 



(36) 



In particular cases (i) and (?z) considered above, Eq. (p9) reduces to expressions 



T, 



cO 



4 T e 



(37) 



and 



8r„ 



(38) 



for initial T c suppression by magnetic (at x — 0) or nonmagnetic (at arbitrary value of x) scatterers respectively. 

As to the high impurity concentration, we recall that in the BCS theory, nonmagnetic scattering alone is insufficient 
for the non-d-wave two-dimensional superconductivity (0 < % < 1) to be destroyed completely ]2(|; at = 0, the 

value of T c asymptotically goes to zero as l/r„ increases. On the other hand, T c of a d-wave superconductor with 
X = 1 vanishes at a critical value l/r„, c = 7rT c o/7 ~ 1.764T c o, with 7 = e c ks 1.781, where C is the Euler constant. In 
its turn, magnetic scattering in the absence of nonmagnetic scattering (1/t„ = 0) is known to suppress the isotropic 
s-wave superconductivity with x = at a critical value l/r^ c = ttTco/2^ w 0.882T c o (Ref. [^7|). 

On the basis of Eq. (|33|), it is straightforward to derive the general condition for impurity (defect) suppression of 
T c for a superconductor having an arbitrary in-plane anisotropy coefficient x an d containing both nonmagnetic and 
magnetic scatterers: 



T eff,c 7 

where T e ff yC is the critical value of the effective relaxation time T e ff, defined as 



1 



1 



i-x 



1 



1 

pot 



1 



(39) 



(40) 



From Eqs. ( p9[ ) and ( [40| ) one can see that l/r e // iC increases monotonically with l/r„, l/r^ ', and 1/r^f at any 
value of Xi with the exception of the case x = 0, where l/r e ff iC doesn't depend on 1/t„ and I/t^ ', see Eq. (f4(i|). If 
X is close to unity (A(p) with strong in-plane anisotropy), then l/r e // « l/r„ + I/t^ ' + , i.e., the contribution 
of nonmagnetic and magnetic scattering to pair breaking is about the same. If x << 1 (almost isotropic A(p)), then 
l/r e ff « 1/r^f, i.e., r e ff is determined primarily by magnetic scattering. The higher the anisotropy coefficient x, the 
greater is the relative contribution of nonmagnetic scatterers to T c suppression as compared to magnetic scatterers. 

We note however that while the concept of the effective relaxation time r e ff can be used for evaluation of the 
critical level of nonmagnetic and magnetic disorder, it is not possible to express T c in terms of r e ff in the whole range 
< T c < T c o, see Eq. ([33]). In other words, the combined effect of nonmagnetic and magnetic scattering on T c cannot 
be described by a single universal parameter depending on the values of t„ , T^ ot , r^f , and x> see Ref. [M for more 
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details. Hence, while the quantity l/r e //,c characterizes the critical strength of impurity scattering corresponding to 
T c = 0, the quantity X/T e tf (when it is less than l/r e // jC ) doesn't determine the value of T c unequivocally. 

Based on Eqs. (^) and (Eoj), it is possible to derive the following expression for the critical value of l/r„ in the 
presence of magnetic scattering: 

1 

This expression is valid as long as its right-hand side is positive, since otherwise the superconductivity is completely 
suppressed solely by magnetic impurities. The value of l/r„ )C decreases as l/rf" 1 and increase at constant \ or 

as x increases at constant l/r^f* and 

To conclude this Section, it is interesting to note that T c doesn't depend on x provided that = 1/t„ + l/r^f*, 

seeEq. (§|). 



/ -kT c0 t% 
V 2 7 



1 

pot ' 



(41) 



4. DISCUSSION 



Equation (|33|) is obviously more general than Eqs. ( p4|) and (|35|), which are commonly used for the analysis of 
experimental data on T c suppression by defects and impurities in HTSCs, see references in the Introduction. In fact, 
making use of Eq. ( |34| ) or Eq. ([35]) one assumes a priori that either (i) the order parameter in HTSCs is isotropic 
in momentum space, or (ii) magnetic scatterers are completely absent in HTSCs. In our opinion, the experimental 
dependencies of T c versus impurity concentration or radiation dose should be analyzed within the framework of the 
theory presented above, see Eq. (|33|). One should not guess as to the degree of in-plane anisotropy of A(p) and the 
type of scatterers, but try to determine the value of x an d relative weights of magnetic and nonmagnetic components 
in electron scattering through comparison of theoretical predictions with available or specially performed experiments. 

We recall that Eq. ( |33| ) has been derived within the weak-coupling limit of the BCS model. Note however that the 
exact solution of the Eliashberg equations for a particular case of a <i-wave superconductor containing nonmagnetic 
impurities only indicates pc| that the analytical T c /T c q versus l/r„ curve falls near the numerically calculated T c /T c q 
versus 1/t* curve, where 1/r* is the scattering rate renormalized by the strong-coupling effects (it is 1/r* that 
enters the formula for the experimentally determined in-plane residual electrical resistivity po). We believe therefore 
that Eq. ([33]) is also valid beyond the weak-coupling approximation implying that r n , r^ *, and r^f in Eq. ( |33| ) 
are the renormalized relaxation times which govern the experimentally measured physical quantities. It would be 
interesting to check this by direct numerical solution of the Eliashberg equations for an anisotropic superconductor 
with nonmagnetic and magnetic impurities. 

In order to compare the predictions of theory with experiment, it is convenient to represent the electron scattering 



time, Eq. (24), in terms of the in-plane residual resistivity pq. Following Radtke et al. [p(|, we have 



1 1 1 u li 

- + ^i + ^ = t P0 > (42) 

In T rn l m 1" 

where oj p i is the plasma frequency. Note that spin-independent (1/t„ + 1/t^ ') and spin-dependent (1/r^f ) scattering 
rates variously appear in Eqs. ( p3| ) and ( [42| ) for the critical temperature and residual resistivity. Hence, for a given 
degree of anisotropy of the order parameter (i.e., for a given value of x), the universal dependence of T c /T c q on po 
cannot be obtained, as opposed to the case of a c?-wave or anisotropic s-wave superconductor containing nonmagnetic 
impurities only [|60| , [3l"l| . 
Let us express po as 

Po = Pr+PT, (43) 

where /?q™ is due to electron scattering by nonmagnetic impurities and potential scattering by magnetic impurities, 
while Pq x is due to exchange scattering by magnetic impurities: 

f + ^ = ^o"\ (44) 

T„ Tfr, 47T 



(45) 
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From Eqs. (E5 



t) we have 



nm 

Po 



a)Po- 



ex 

Po 



ap , 



where 



,ex\2 



a 



(c n /c m )\u n \ 2 + \U 



pot, 2 



,ex\2 



(46) 



(47) 



The value of a depends, first, on the scattering strengths of individual nonmagnetic and magnetic impurities (through 



matrix elements 



pot 



c ) and, second, on the ratio of impurity concentrations c n /c m . The latter is expected to 



remain constant under doping or irradiation, at least at relatively low (but sufficient to destroy the superconductivity) 
doping level or radiation dose. For example, low energy irradiation of YBa2Cu307_a; was found to induce nonmagnetic 
defects only [ j50[ , i.e., c m /c n = 0, and hence a = 0. 

Thus the dependence of T c /T c q on po for a given value of x is specified by the material-dependent and " disorder- 
dependent" dimensionless coefficient a. The greater is the relative contribution from exchange scattering by magnetic 
impurities to po, the higher is the value of a (a ranges from in the absence of exchange scattering to 1 in the absence 
of non-spin-flip scattering). Substituting Eqs. (Ej) and (W3) in Eq. (B3) and taking Eqs. (Efl) into account, we have 



In 



rO 



= (i-x) 



1 



8tt 2 T ( 



Po 



- ^ 




16tt 2 T ( 



Po 



- * 



(48) 



Figures 1-4 show the plot of T c /T c q versus po in a superconductor with T c q — 100 K and ui p i = 1 eV for different 
values of \ and a ranging from to 1. The choice of T c q and oj p i is, to some extent, arbitrary (though these values 
of T c o and ui p i are typical for HTSCs, e.g., for YBa2Cu307). In order to go to the other values of T c q and oj p i one 
should just replace po in Figs. 1 - 4 by po(T c o/ 100)u;~ ; 2 , where T c q is measured in K, and uj p i is measured in eV. 

From Figs. 1-4 one can see that at % < 1 the rate of T c decrease with increase in po becomes higher as a increases 
from to 1, i.e., as the relative contribution of exchange scattering to po increases. At \ — (isotropic s-wave pairing) 
the value of T c does not depend on po for a — 0, while the superconductivity is completely suppressed (T c = 0) at 
a critical value of pg = 1.42 mficm, 113 p,£lcm, and 56.5 pSlcm for a = 0.04, 0.5, and 1 respectively, see Fig. 1. At 
X = 0.5 (a specific case of anisotropic s-wave or mixed (d+s)-wave in-plane pairing) the value of T c monotonously goes 
to zero as po increases for a — 0, while p[j — 401 pQcm, 113 pQcm, and 80 pflcm for a — 0.04, 0.5, and 1 respectively, 
see Fig. 2. At \ — 0.8 (strongly anisotropic s-wave or mixed (d + s)-wave in-plane pairing with predominance of 
d-wave component) one has p^ = 188 pftcm, 113 pilcm, and 99 pVtcm for a = 0.04, 0.5, and 1 respectively, see Fig. 
3. The curves T c (po) for different a come closer together as the coefficient \ increases, i.e., as the order parameter 
becomes more anisotropic. At % = 1 (d-wave in-plane pairing) all curves T c (pg) merge together, see Fig. 4, i.e., the 
value of T c /T c o at a given po does not depend on a, in accordance with Eq. (|48|), the critical value of p§ being equal 
to 113 pilcm at any a. Note that for a = 0.5 the curves T c (po) are the same at any value of x, see Figs. 1-4 and 
Eq. @. 

Magnetic scatterers in a non-d-wave superconductor, even if they are present in a small proportion (a << 1), result 
in p^ decrease as compared with of a sample containing nonmagnetic impurities only. The decrease in p^ with a 
is more pronounced at low values of Xi i- e -> m superconductors having weakly anisotropic order parameter, see Figs. 
1-3. At % as high as 0.8, i.e., in a superconductor having strongly anisotropic (but different from a pure d-wave) 
order parameter, the value of for a = 1 is less than twice as low as that for a = 0.04, see Fig. 3. In such a 
superconductor, the role of a small amount of magnetic impurities is to suppress the superconductivity at a finite 
value of p§ as opposed to the case when exchange scattering is absent (a = 0), though the curves T c (po) at a — and 
0.04 almost coincide in a very broad range of T c /T c o, see Fig. 3. 

In our opinion, an argument in favor of other than pure d-wave in-plane symmetry of the order parameter in HTSCs 
(at least in some of them) is as follows. A pure d-wave two-dimensional superconductor with x = 1 is characterized 
by the universal dependence of T c on po which is the same at any value of a, i.e., at any relative contribution of 
exchange scattering to the total value of po, see Eq. (^8|) and Fig. 4. Meanwhile, T c versus po curves and the values 
of p% in HTSCs are material-dependent and disorder-dependent Q||,||j44|j4|j5|,||,|5|j5^-|9| . This fact attests that 
the value of x varies (though, may be, slightly) from one HTSC to another, while the value of a depends both on the 
kind of HTSC material and on the type of impurities or radiation-induced defects. 

Besides, the experimentally observed form of T c (po) curve in HTSCs is usually close to linear in a very broad range 
of critical temperatures (42|4|^,||,||,|5|,||-|9| . The theoretical curve T c (po) has such a form if x is close to unity 
(but x 1) an d ck is much less than unity, e.g., at x = 0-8 and a = 0.04, see Fig. 3. In contrast, the theory predicts 
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the negative curvature of T c (po) curve for a pure d-wave superconductor (no matter how great is the contribution of 
exchange scattering to po), see Fig. 4, as well as for a non-d-wave superconductor with strong exchange scattering, and 
the positive curvature of T c (po) curve for a non-<i-wave superconductor containing nonmagnetic impurities only, see 
Figs. 1-3. So, we expect that the majority of HTSCs have the mixed (d+s)-wave order parameter with predominance 
of d-wave component (1 — % << 1) and that exchange scattering by magnetic impurities or radiation- induced defects 
contributes to po, though quite insignificantly (a << 1). 

The admixture of s-wave component to a <i-wave order parameter, e.g., A(p) = Ad{cosp x — cosp y ) + A s , may 
be a consequence of orthorhombic distortion of CuC-2 planes in some HTSCs jl7]-jl9|. The value of the coefficient x 
contains the information about the partial weight of that component in the order parameter, i.e., about the value 
of A s /Ad. So, having determined the value of x from experimental data on radiation-induced and impurity-induced 
reduction of the critical temperature, one can deduce the value of A s /Ad making use of Eq. (pl"|). 

Besides, it should be stressed that (d + s)-wave symmetry is only one of possible candidates for the symmetry of 
anisotropic pairing state in HTSCs. It is likely to occur in orthorhombic HTSCs. In what concerns purely tetragonal 
HTSCs, one may expect mixing of isotropic s-wave state with the state having some higher even angular harmonic, 
e.g., with g- wave state. Such a mixed {g + s)-wave state, just as (d + s)-wave state, is also characterized by x values 
in the range from to 1, depending on the partial weights of s-wave and g-wave components in the order parameter. 
All the results obtained in this paper are therefore applicable to the case of (g + s)-pairing, as well as to the case of 
any other in-plane symmetry of the order parameter. 

To conclude this Section, we note that an assumption about the constancy of the parameter a (i.e., an assumption 
about the constancy of the ratio of the concentrations of nonmagnetic to magnetic scatterers) under doping or 
irradiation must be checked before detailed comparison of the theory presented in this paper to experimental data. 
If this assumption appears to be incorrect, Eq. (|33|) for the critical temperature can still be used, the scattering 



times being given by Eqs. (44) and ([45|). In that case, however, one faces an additional complication concerning 
the evaluation of contributions to the residual resistivity po from magnetic and nonmagnetic scatterers, Pq x and p$ m 
respectively. 

5. SUMMARY 

The combined effect of nonmagnetic and magnetic defects and impurities on the critical temperatures of supercon- 
ductors with different gap anisotropies was studied theoretically within the weak coupling limit of the BCS model. 
For the case of short-range scattering potentials, an expression was derived which relates the critical temperature to 
the relaxation rates of charge carriers on nonmagnetic and magnetic scatterers as well as to the coefficient of in-plane 
anisotropy of the superconducting order parameter on the Fermi surface. 

We note that the results obtained in this paper can be modified to include the effects of anisotropic (momentum- 
dependent) impurity scattering. For example, in the case of significant overlap between the anisotropy functions of 
scattering potential and that of the pair potential, the anisotropic superconductivity has been proven to become less 
sensitive to nonmagnetic impurities [|62| , |63"|]98f| . However it is not clear if there is such an overlap in HTSCs. 

Besides, numerical calculations within an extended Hubbard model point to the spatial variation of the order 



parameter in the vicinity of impurities in anisotropic superconductors |66|. As a result, suppression of T c is significantly 
weaker than that predicted by the Abrikosov-Gor'kov-type theory. This effect presumably is especially pronounced in 
superconductors with short coherence length. However, a complete theory of such an effect remains to be developed. 

It is worth noting that impurity doping and irradiation generally result not only in a structural disorder but also 
in creation or annihilation of charge carriers. Thus the effects of carrier and impurity concentrations on T c of HTSCs 
should be considered on equal footing |99[] . Moreover, since high-temperature superconductivity appears upon doping 
of parent insulators, a description of those effects should form the basis for the future theory of HTSCs. 

In conclusion, the results obtained provide a basis for evaluation of the degree of anisotropy of the superconducting 
order parameter (e.g., for an estimate of the partial weight of s-wave in mixed (d+ s)-wave order parameter) as well as 
the ratio between nonmagnetic and magnetic scattering rates in high-T c superconductors through careful comparison 
of theoretical predictions with the experiments on impurity-induced and radiation-induced reduction of the critical 
temperature. We hope that the present paper will serve as a stimulus for further experiments on combined effect of 
nonmagnetic and magnetic scattering in the copper-oxide superconductors. 
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FIGURE CAPTIONS 



Fig. 1. Dependence of the normalized critical temperature T c /T c q on the residual resistivity po due to nonmagnetic 
and magnetic impurities in a superconductor with T c q = 100 K and x = (isotropic s-wave pairing) for different 
values of the coefficient a specifying the relative contribution to po from exchange scattering, a — (solid curve), 0.04 
(long-dashed curve), 0.5 (short-dashed curve); 1 (dot-dashed curve). The critical value p% — 1.42 mVlcm for a — 0.04. 
The plasma frequency is taken to be tu p i = 1 eV. One can go to the other values of T c q and lo p i through replacing po 
by p (T'co/100)(jj~ ; 2 , where T c0 is measured in K, and uj p i is measured in eV. 

Fig. 2. Same as in Fig. 1 for \ = 0.5 (a specific case of anisotropic s-wave or (d + s)-wave in-plane pairing). 

Fig. 3. Same as in Fig. 1 for x = 0.8 (a specific case of anisotropic s-wave or (d + s)-wave in-plane pairing). 

Fig. 4. Same as in Fig. 1 for x = 1 (d-wave in-plane pairing). In this case the value of T c /T c q at a given po does 
not depend on a, see Eq. (fi"8|). 
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